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Low frequency shot noise and do current profiles for a double-barrier 
resonant-tunneling structure (DBRTS) under a strong magnetic field applied 
perpendicular to the interfaces have been studied. Both the structures with 
3D and 2D emitter have been considered. The calculations, carried out with 
the Keldysh Green's function technique, show strong dependencies of both the 
current and noise profiles on the bias voltage and magnetic field. The noise 
spectrum appears sensitive to charge accumulation due to barriere capaci- 
tances and both noise and dc-current are extremely sensitive to the Landau 
levels' broadening in the emitter electrode and can be used as a powerful tool 
■ to investigate the latter. As an example, two specific shapes of the levels' 

broadening have been considered - a semi-elliptic profile resulting from self- 
consistent Born approximation, and a Gaussian one resulting from the lowest 
0^ ' order cumulant expansion. 

(N 

(N 
> 

^ . I. INTRODUCTION 

' In recent years, there has been a great interest in resonant tunneling through double- 

I barrier resonant tunneling structures (DBRTS) (Fig. 1). Such structures have been in focus 

of many experimental and theoretical investigations since its conceptioruby Tsu and Esaki^l 
and first realization of negative differential resistance by SoUner et aln. Many important 
characteristics of DBRTS have been intensely studied, e.g. dc-properties , phonon assisted 
tunneling, time dependent processes and frequency responsej-. Noise properties of DBRTS 
have also been studied both experimentallyB and theoreticallyacl. At low temperatures and 
^ ' in the presence of transport current, shot noise is the dominant source of electrical noise. 

Q , This kind of noise is due to discreteness of the electron charge, and it is sensitive to the 

O ■ degree of correlation between tunneling processes. In general, a correlation leads to an 

additional frequency dependence of shot noise, as well as to its suppression below the so- 
called full noise, S{0) = 2e|/dc| (at T = 0)0. Here S{uj) is the noise spectrum (see the exact 
definition below), while Idc is the average dc current. In a mesoscopic conductor having 
^ . several independent modes of transverse-mjotion (channels) , the noise is determined by the 

partial transmission probabilities astEltj S oc X)m^"i(l ~ Tm), while the conductance 
goes as G oc J^m ■ Suppression of the shot noise is thus expected in a phase coherent 
system when the tunneling probabilities are of the order unity for open quantum channels. 

Our concern is a DBRTS in a strong magnetic field perpendicular to the interfaces. Mag- 
netic field is an important tool for sample characterization because it leads to the formation 
of Landau levels, as well as to drastic modification of electron wave functions. We study 
the situation when the magnetic field B is applied parallel to the tunneling current I, as 
schematically illustrated in Fig. 1. In such a configuration, the magnetic field leads to an 
effectively one-dimensional tunneling problem. Consequently, both the dc current and the 
noise appear extremely sensitive to the details of the density-of-states behavior. We believe 
that such a sensitivity can provide a powerful tool to investigate details of the Landau levels' 
broadening in resonant tunneling structures. 

The paper is organized as follows: Section || describes the model Hamiltonian as well 
as the basic expression from which the current and shot noise profiles will be derived in 
Section [II. In Appendix ^ and Appendix ^ the Green's functions used in our calculations 
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are expanded. 
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As an example, we consider a GaAs+/ Alo.3Gao.7As/ GaAs/ Alo.3Gao.7As/ GaAs"*" 
DBRTS, with the barriers' and the well widths of the order of 40-60 A. Such structures 
were extensively studied experimentally. In many cases the barrier height is about 300 
meV, and it is assumed that there exists only one quasi-bound state in the well. 



II. MODEL AND BASIC EXPRESSIONS 



Consider a DBRTS in the presence of an external magnetic field B perpendicular to the 
interfaces which are assumed perfect, B || I || z. Within the quantum well, the electron wave 
function can be expressed as a product of a quasi-bound state xi^) times a wave function 
correspondent to the motion in the x — y plane. Let us denote the energy of the motion 
in z-direction as eo. Under the Landau gauge A = (0,i?x, 0) the wave functions can be 
specified by the set of quantum numbers a — (n, ky) as 

0a(r) = -T= exp{ikyy)ipn{x + Pky)x{z). (1) 
vLy 

The corresponding energy levels (measured from the conduction band edge) are 

E^ = E„^eo + huj^{n+^). (2) 

Here, ipn{x) denote harmonic oscillator states, Uc = eB/m* is the cyclotron frequency, and 
I = y^Ti/eB is the Landau magnetic length. 

Similarly, electron states in the leads are specified by the quantum numbers P = 
{m, kj^y, kj_z), where j = e(c) refers to emitter(collector) states, respectively. The corre- 
sponding wave functions and energy levels under the bias eV are given as: 

(l)j,fi{r) ^ exp{ikj^zZ + ikj,yy)ipm{x + l^kj,y), (3) 



where < ag < 1 and Ac = flg — 1 (the symmetric = 0.5 case will be considered in our 
numerical calculations). We arrive at the model Hamiltonian 

H = 'Hc+'Ht^^ Ej,l3c]pCj^l3 + ^ EaC^^Ca + ^ [V,-,/3a4cj,/3 + h.c] , (5) 
3,13 « 3,a,P 

where the tunnel matrix elements Vj^pa have to be calculated using the eigenstates listed 
above. Since the interfaces are assumed to be perfect, the quantum numbers n and ky 
are conserved during the tunneling process, and so the calculation oLthe matrix elements 
Vj.j3a reduces to the solving of a one-dimcpsional Schrodinger equationEj, following with the 
application of the Bardeen's prescriptioriiJ. Consequently, the tunneling matrix elements 
can be written as 



Vj,Pa — Sm.n^{ky kjy)Vjn{kj^y,kj;,)- (6) 

In noise calculations the time dependence of the tunneling currents flowing through the 
DBRTS is important, and hence the junction capacitances should be taken into account. 
The effect of the junction capacitances can be inciujded in our model with the help of an 
equivalent circuit of the DBRTS as shown in Fig. 2.li3 In this circuit, we specify the currents 
through the emitter (collector) barriers /o(c)(0 and their resistances as i?c(c)- The "external" 
current I{t) is in this model given by 

i{t) = ^m + ^m, (7) 
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where Co(c) is the capacitance of the emitter (collector) barrier and C = Cc + Cc is the total 
capacitance of the quantum well. In the symmetric case Ce = Cc, we arriwe at the simple 
relation I{t) = [Ie{t) + Ic{t)]/2, which was the basis of the Chen & Ting'so calculation for 
shot noise in a DBRTS in a,zero magnetic field. If one ignores the charge accumulation, all 
three currents are the sameO, I[t) = Ic{t) — Ic{t) and the result in this case can be obtained 
from the following formulas in the limit of strong asymmetry, Cc(c)/C — » 0. The asymmetry 
in capacitances is of course not important for the dc-current, where 

Idc = Idc,e = Idc.c (8) 

In the further analysis it is convenient to out ?i = 1, and then restore h again in the final 
expressions and order-of-magnitudc estimates. 

The tunneling current flowing into the well from the emitter and the current /c flowing 
out of the well to the collector, are in general different. They are given by the expressions 

= -i2e^, J2 [VjMci{t)cjAt)) - VlpMA*^Ca{t))] (9) 

I3,a 

where Nj{t) are the Heisenberg number-of-particles operators. Kg = 1, Kc = and a spin 
degeneracy factor 2 is introduced. 

The shot noise spectrum is defined as the Fourier transform of the current-current auto- 
correlation function asu 

/OO /•OO 
S{t)e"^'dt = 4 / S{t) cos{ujt)dt (10) 
-OO Jo 

where S(t) is the quantum mechanical and statistical average of the current-current anti- 
commutator: 

Sit) = 1 {{Am , A/(0)}) = 1 ({/(t) , /(O)}) - II. (11) 
From (^ and (|^) , it can be expressed (having in mind the spin degeneracy factor of 2) as 

Sit)^-e^ J2 VjVjo^ [l^J,/3aV,-o,/3oao ({4WCj,/3(t),cl„(0)Cjo,^;„(0)}) 



- yj*f3aVjo,Poao ( { 4/3 W ' ^lo WCjo , A, (0) } ^ 

+ ^:^o^;o,,o.o ({4,(0c.(t),4_,,,(0)c„„(0)})] . (12) 

where rje = Cc/C and 7]c = —Ce/C. Being expressed through Feynman's graphs, these 
averages involve only the diagrams with the Green's functions connecting the times t and 0, 
since disconnected parts are all canceled by the /J^, subtraction in (|ll|). 



III. THE RESULTS 

The task is now to expand the quantum statistical averages appearing in (^) and (p^. 
For a finite bias, the DBRTS as a whole is not in thermal equilibrium, and it. SRems thus 
appropriate to employ the Keldysh non-equilibrium Green's function techniqueEjtJ, where 
the two lead subsystems are supposed to have their own local equilibrium. 

Expanding (^ yields (Appendix]^) 

/•OO 

^dcj = E / lj{n,e)A{n,e) [fQw{e) - fj{e)] . (13) 
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In the above expression, gi, = Ly/2TTp is the magnetic ky summation degeneracy factor, 
jj{n,e) is the escape rate to the lead j, fgwis) and fj{e) are the occupation factors, while 
A{n, e) is the spectral function for nth Landau level in the well, 

Ain,e) = -23Gfi(n,e) = ^ (14) 

Here, GR{n,e) is the retarded electron Green's function, and 7(n,e) = 7e(ri,e) +7c(n,e) is 
the level broadening due to the finite escape rate to the leads. Usually, the energy distance 
between the resonant level in the well and the tops of the barriers is much greater than 
the escape rate from the well, 7. In this case the tunneling matrix element (|^) can be 
considered as a smooth function of the energy in comparison with the energy dependence of 
the density-of-states in the leads, 
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Thus the escape rates 7^ can be expressed as 

7j(n,e) = 2Tr\Vj\'^gj{n,e). 

Consequently, the noise appears a sensitive tool to study density-of-states in the electrodes 
in a magnetic field. Below, we will do numerical calculations for two models for the density 
of states - for a constant Lorentzian broadening, and for the so-called self-consistent Born 
approximation. 

Since both leads are assumed to be in a thermal equilibrium with different electro-chemical 
potentials and Fermi energies, the occupation numbers can be expressed as the Fermi func- 
tions: 

^^^^^ " ^ie-Ef-a,eV)/kBT ^ ^ (1^) 

However, thermal equilibrium is not maintained in the quantum well and thus one cannot 
use the Fermi distribution for the electrons in this region. Instead, from tte dc-current 
conservation law (^, the weighed average occupation factor is determined asliSl 

f I X 7e(n,£)/e(£) +7e(n,£)/e(£) , , 

fQwin, e) = -^^^ . (16) 



Rc-introducing h to return to the proper units, we arrive at the Landauer formuk 
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/dc = ^ 5] / d£ T„(£) [Ue) - Ue)\ (17) 
with the transmission probability 

T„(£)^ "-("'f-^;'^) A(n,£) (18) 
7(n,£) 

To get a relatively simple expressions for the shot noise from Eq. (|l|) we assume the 
following approximations. First, we assume that the resonant level is situated well inside 
the resonant tunneling region, 

\ajeVj + Ep — En\ ^ max(?i(jj, 7, Vj), 

\ajeVj ~ Col > max(?itj, 7, z^j), (19) 

and that lo <C Wc- These inequalities allow us to put /j(£ ± w) — > jj{s) and 7j(n, £ ± w) — > 
7j(n, £). Second, the temperature is assumed to be low {ksT ^ 7), in which case the Fermi 
functions can be approximated as step functions. Keeping those approximations in mind, 
we arrive at the following result (Appendix P) : 
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"7e7c(7eCe - 7cCc)(7eCc - 7cCe) 7e7c"' 



X e)A{n, e — uj) 

+ [A{n,e) + A{n,e- uj 



C272 72 

C2 7 



- m [Gnin, £)] [Gnin, e - w)] ^^jelJj (20) 

where 7^ = 7j(n,e). 

Re-inserting h, and using the relation Ajelc^ [GR{n, e)]^ — 4T„(£) — (7^/70 7c) T^i^), we 
arrive at the well known result: 

^(0) = ^ / de T„(e) [1 - T„(e)] [/e(e) - Ms)]' (21) 

As one could expect, the zero frequency shot noise does thus not depend on the barriere 
capacitances and the aboyejE^ult coincides with previous calculations which have been, per- 
formed for point contactsElEJ, for arbitrary phase coherent two terminal conductorsEj (ne- 
glecting barriere capacitances) , and also for a DBRTS in the regime of incoherent tunnelingp. 
The main features of our problem is that the combinations T„(l — T!„) enter for each Landau 
level independently and that the tunneling probabilities T„ are strong functions of magnetic 
field. An important feature is that Eq. (|T[) holds even if the inequality ( p^ is violated. That 
makes zero-frequency shot noise, together with the dc-current, a powerful tool to investigate 
the density of states in the leads which manifests itself through the escape rates jj . 

The results for a particular DBRTS device are shown in Fig. 3. Here we use the model 
of constant Lorentzian broadening of the Landau levels, where the escape rates can be 
expressed as (Appendix ^ 

7j(n,£) = — 



2^/2 {e- Ej,n)' + {i^/2)' 



Here, T^- is a constant characterizing the strength of the escape rate and = eV (ij + 
uic^n -f 1/2). Note that there are peaks in the dimensionless shot noise factor S{ijj)/eldc at 
the voltages when an intra-well Landau level passes the emitter's electro-chemical potential. 
Those peak's shape is determined by an interplay between the quantum suppression (S'(O) oc 
T„(l — T„)) and a finite broadening of the Landau levels in the quantum well. In addition, 
a small peak appears in the dc-current curve at the end of the resonant tunneling region (in 
our example, at eV ~ 55 meV) due to the finite broadening of the lead electron states. This 
broadening can typically be of the size v ^ he/m* fi ~ 0.5 meV (/i is the electron mobility). 

The effect of the level broadening in the leads is even more pronounced in the case of a 
2D-emitter. For numerical calculations in this case we employ the so-called self-consistent 
Born approximationE3. In this approximation, the density of states takes a semi-elliptic form 
and the escape rate is then given by: 



where E^.n = eVae+LiJc{n- + l/2) + ee, Se is the emitter quasi-bound level and Lez is the width 
of the 2D-emitter. The lead broadening depends in this case on the magnetic field and is 



given hy V ^ \j2Ti e ujc/i^m*^, where /z is the mobility of the 2DEG. In our example, ~ 10 

cm^/Vs, at huoc — 10 meV we get v ^ 0.35 meV. In realistic systems, sharp edges of the semi- 
elliptical densitaof-states profile are smoothed, the smoothing for a long-range potential 
being GaussianEil. To check the sensitivity to the smoothing we made also calculations for 
a Gaussian density-of-states profile. The calculations show that both the current and the 
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noise profiles can be very sensitive to the degree of such a smoothing. Fig. 4 shows the dc- 
current and zero frequency shot noise results for a particular DBRTS device with 2D emitter 
calculated according to the self-consistent Born approximation (semi-elliptic profile) as 
as for a Gaussian profile obtained from a so-called lowest-order cumulant approximatior 
A double-peak structure is obtained with the Gaussian profile in contrast to the single peak 
appearing in the case of a semi-elliptic profile. 

We believe that our results can serve as a basis for an experimental test of the strength 
of the Landau level's smearing by impurities. In our example, the splitting of the noise and 
current peaks in the case of the Gaussian level broadening case is about 2 meV, and should 
be observable at temperatures T ^ 20 K. 

Finally, we give an expression for the shot noise valid at finite frequency provided the 
inequality (19|) holds. Integrating (|l^ and (^) with a 3D emitter, we arrive at the expression 
(Appendix 



Siu;) = 



2\eldc 

fJ2 



1 



[CeC,(7e'+7c)-(C^e+Ce')7e7c 
7e7cC2 -t- CeC,7'] } 



(24) 



This result is strongly dependent of the bias voltage because of the voltage dependence of 
the escape rates. Indeed, at |eo — eyaj j 3> max(?ia;, 7, j/j). 



1] = 7j(eV") 



8(eo — eVaj 
yJ(.Q - eVaj 



As our two special cases, symmetric capacitance (Ce = Cc) and no charge accu 
(Cc(e) ~f-|0), we arrive at the relations similar to those obtained by Chen & Ting 
BiittikerO in zero magnetic field: 



eld 



7^ -I- a;2 



,7e7c 



dc 



1 



7e7c 



7^ -|- uj^ 



iiulation 
i and by 



(25) 
(26) 



However, the important difference is strong dependencies of the escape rates on both electric 
and magnetic fields. The frequency dependency of the noise in those two cases are very 
different (Fig. 5) and can serve as a basis for an experimental test of the importance of the 
charge accumulation on the barrier capacitances in the DBRTS tunneling structure. 

The present work has partially been supported by the Norwegian Research Council, Grant 
No. 100267/410. 



APPENDIX A: GREEN'S FUNCTION EXPANSION FOR DC CURRENT 

The quantum statistical averages appeaxing in (^ is expanded using the Keldysh non- 
equilibrium Green's function techniquelljO. Four different Green's functions, appropriate 
for a S-matrix expansion in the time-loop formalism, are defined along a closed time path 
that runs from — oo to -l-oo along the a —'V branch and then returns from -foo back to — oo 
along the a — '2' branch: 

G„ia. ih - h) - -z(7^c(ti)ct(t2)) (Al) 

where by (t„ = 1(2) is meant that i„ is located on the T'('2') branch and % is the generalized 
chronological operator ordering physical operators along the closed time path. In the Fourier 
transformed energy space, the Green's functions are simply related to the retarded Green's 
functions as: 

Gn{s)^if{e)A{e) + GR{e), 
Gi2(e)-z/(e)A(e), 
G2i(e) = -z [l-fie)]Aie), 

G22{e) = -i[l- /(e)] A{e) - GR{e). (A2) 
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Here A{e) = —29 [G'fl;(e)] is the spectral function, while f{e) is the occupation number in 
the region considered. The following retarded quantum well and lead Green's functions are 
used as the basis in the calculations: 

G%{jP, E,jj + iv,l2Y\ (A3) 

Here 7(n,£) = 7e(r7,,e) + 7c(n,e) is the broadening of the resonant states due to the finite 
tunneling rate to the leads, and Vj is the broadening of electron states in the leads due to 
electron scattering. |— , 

The dc current is expanded to lowest order in the time-loop S-matrix expansionllZI, which 
from (^) yields (as diagrammatically represented in Fig. 6): 

/dcj = -4eKj lim V / dti I V,,0^ |2 SR IrA Ati)co,{tx)c,^fi{t)c^^(i) 

/OC 



~Gi2{jl3,t-h)G22{a,h-t)]. (A4) 



In the first of the above integrals, t is located on the '1' branch, t is located on the '2' branch 
and the ti integral is taken along the time-loop from —o o to -t-oo and back to — cxd. The 



latter result, introducing Green's functions according to (Al), is expressed as an integral 



over ordinary re al t ime axis from — oo to oo. The Fourier transform of this result, with the 



substitution of (A2), yields: 



/oo 
de A{a,e)A{j(3,e) [fQw{e) - /,(e)] , (A5) 

where fqwi^) f^nd fj{e) are respectively the occupation numbers in the quantum well and 
leads. Using the escape rates from the quantum well states to the lead j, defined as 



7j 



(n,s)=^|F, |2A(j/?,£) (A6) 



where the tunneling matrix elements in (^ have been assumed to be independent on any 
quantum numbers (Vj,n(fcy,fci) = Vj) and taking into account the ky independence of the 
electron Green's functions, {A{a,s) = A{n,e)), we arrive at ( p^ ) and (p2|). 



APPENDIX B: GREEN'S FUNCTION EXPANSION FOR SHOT NOISE 

The quantum statistical averages appearing in (^2|) is expanded in a similar way as with 
the dc current. It is found that S{uj) is symmetric in uj and can be written as a sum of 6 
different terms (represented by the diagrams in Fig. 7): 

S{u;) = Si-tj) = Siiu;) + 52(c^) + 83(0;) + Si{ij) + S^itj) + Seicu). (Bl) 

Si{uj) is expanded from the first term in ( p^ ) as 

Si{uj) ^ Siai^) + Siai-^) (B2) 

with 

=2 



j,jo,a.0 rTi,cr2 

X [Ga2i{a,e)G2cr2{jf3,s)Gai2{a,e - uj)Gicri{jol3,£ - i^)] ■ (B3) 
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S'2(w) is the contribution from the from the 4'th term in ([T^), simply related to Si{lu) as 

52 (w) = Sliuj). (B4) 
The second term in (f2|) has both zeroth [53(0;)] and second order [5*4(0;)] contributions: 



53(0;) 
54(0;) 



e 

TT 



e 

TT 



r,^ \Vj \^ J deG2iijf3,e)Gi2ia,e-Lu) + Gi2UP,e)G2i{a,s-u;) 



CToCTt = {l2,2l} <Tl,<T2 

X [^,71^70 £)Ga2a, (a, e)G'<TtfT2 (jo/3, e)^,^^^, {a, e - uj)] . 



+02 



(B5) 



The zeroth [55 (tj)] and second order [56 (w)] contributions from the third term in ( |12D are 
simply given as: 



55(w) = S-i{-uj) 
Sq{uj) = Si{-uj). 



(B6) 



The diagrams of the type shown in Fig. 8 are not taken into account explicitly because 
they are already included in 53 (w) and 55(0;), since the quantum well electron Green's 
functions we use as our basis are originally dressed by tunneling to the leadsQ. Summing 
up the different diagrammatic terms we neglect the contribution from real part of the lead 
retarded Green's functions. This is a reasonable approximation since it corresponds to a 
Hilbert transform of the imaginary part (proportional to the escape rates) and it appears 
that, in and above the resonant tunneling region, its contribution is negligible. Keeping 
this in mind, as well as the approximations listed in the main text (/j(e ± tj) — > /j(e), 
7j(n, e ±uj) ^ Ijin-, e) and ksT ^7), we arrive at the result ([20|). 

Integrating the shot noise expression ( po|) and the dc-current (|l^), we make use of the 
following integrals over the resonant tunneling region, valid for a Landau level located well 
inside the resonant tunnehng region according to (19): 



J deA{n,e) « 2tt 



deA{n, e)A{n, e — lo) 
deRe [Gii{n, e)] Re [GR{n, e — u)] 
With those relations, we arrive at (p4). 



47r7 

72 -j- a;2 

777 
72 + a;2 



(B7) 
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FIG. 1. Schematic illustration of the double-barrier resonant tunneling structure (DBRTS). 
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FIG. 2. Equivalent circuit for a DBRTS. 
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FIG. 3. Average current and zero frequency shot noise (A) and dimensionless noise-to-current 
ratio (B) for a symmetric SD-emitter DBRTS with: Te = = 0.67 (meV)^/^ hujc = 10 meV, 
eo = 27 meV, Ef = 30 meV and v = 0.5 meV. The shot noise ratio sohd curve was obtained from 
the exact equation (pi]), while the dashed curve results from the approximation (H). 
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FIG. 5. Frequency dependence of the noise-to-current ratio for a symmetrical DBRTS with the 
parameters Te = Tc = 0.67 meV'^''^, huic = 10 meV, £c = 27 meV, Ep = 30 meV and with the 
apphed voltage eV = 30 meV. Curve '1' shows the case of symmetric barrier capacitances (|2^), 
while curve '2' is the result when barrier capacitance charge accumulation is negligible (|2^). 
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FIG. 8. A typical diagram not taken explicitly into account since it is already implicitly included 
in other diagrams. 
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